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are  presented  and  the  best  and  worst  fiber  orientations  are  identified 
Finally,  the  interlaminar  stresses  are  computed  and  the  variation 
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Part  II 

This  part  deals  with  the  3D  stress  field  ot  a  cylindrical  fiber  which  is 
embedded  into  a  resin  matrix.  The  composite  is  then  subjected  to  a  uniform 
tensile  load  c0  The  strain  energy  release  rate  is  computed  and  the  criterion  is 
used  to  predict  debonding  initiation  at  the  fiber/matrix  mtertace  The 
analysis  shows  that  this  failure  is  most  likelv  to  occur  at  the  free  surface,  le  the 
region  where  the  fiber  intersects  a  free  surface  for  example  a  hole,  an  edge,  or 
a  crack.  Moreover,  it  will  occur  at  approximately  (1/10)  the  load  value 
required  for  the  same  failure  to  commense  at  the  center  of  the  fiber  length 

The  results  are  also  extended  to  include  a  doubly  periodic  array  of  fibers 
which  are  embedded  into  a  matrix.  Based  on  3D  considerations,  the  stiffness 
matrix  is  shown  to  increase  as  the  volume  fraction  of  the  fibers  increases. 
Similarly,  the  stress  Orr  in  the  matrix  is  shown  to  decrease  as  the  volume 
traction  of  the  fibers  increases. 
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PREFACE 


It  was  found  convenient  for  development  and  clarity  purposes  to  divide 
the  contents  of  this  report  in  two  parts.  Part  I  deals  with  the  3D  stresses 
in  a  laminated  composite  plate  weakened  by  a  circular  hole.  The  analysis 
is  based  on  3D  macromechanical  considerations.  Part  II  deals  with  the 
debonding  aspects  between  a  fiber/matrix  interface  particularly  in  the  re¬ 
gion  when '  the  fiber  intersects  a  free  edge.  e.g.  the  surface  of  a  hole.  Tie- 
latter  results,  which  are  based  on  micromechanical  considerations,  are  then 
used  in  conjunction  with  the  results  of  part  I  to  predict  the  critical  applied 
load  stress  which  may  cause  initiation  of  ply  delaminatiou. 


Part  I 


THE  3D  STRESS  FIELD  ADJACENT 
TO  A  HOLE  L\  A  LAMINATED  COMPOSITE  PLATE 
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Abstract, 


An  eigenfunction  expansion  is  developed  for  the  determination  of 
the  exact  31)  stress  field  in  the  neighborhood  of  the  intersect  ion  of 
tin'  free  edge  of  a  hole  and  an  interface  in  a  laminated  composite 
plate.  For  transversely  isotropic  laminae,  the  st ress  field  is  shown 
to  possess  a  weak  singularity  whose  strength  depends  on  tin  material 
constants,  the  fiber  orientations  of  the  two  adjacent  laminae  as  well  as 
the  polar  angle  0.  Results  for  [0°/y0°] ,  [U  ' / “ O  ’] ,  [O  ’/ -10 “] .  and  [U  /'JO  ’ 
are  presented  and  the  best  and  worst  filmr  orientations  are  identified . 

Finally,  the  interlaminar  stresses  tire  computed  and  tin  variation 
as  a  function  of  the  angle  0  is  identified.  The  circumferential  stress 
(Too  is  shown  to  posses  a  small  jump  across  the  interface.  Results  for 
its  behavior  in  the  interior  of  each  layer  are  also  given. 
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1  Int  rod  action 


Despite  careful  design,  practically  every  structure  contains  stress  concentra¬ 
tion.'  clue  to  holes.  Bolt,  holes  and  rivet  holes  are  necessary  components  lor 
structural  joints.  It  is  not  surprising,  therefore,  that  the  majority  of  service 
cracks  nucleate  in  the  visinity  of  a  hole.  While  the  subject  of  stress  con¬ 
centrations  is  certainly  familiar  to  engineers,  the  situation  is  significantly 
more  complex  in  the  case  of  high  performance  laminated  composite  mate¬ 
rials.  The  presence  of  a  hole  in  the  laminate  introduces  significant  stress 
contributions  in  the  third  dimension  which  create  a  very  complicated  hi) 
stress  feild  in  the  vicinity  of  tin'  hole.  Moreover,  this  complex  state  of  sties' 
may  depend  on  the  stacking  sequence  of  the  laminate,  the  fiber  orientation 
of  each  lamina  a.s  well  as  the  material  properties  of  the'  fiber  and  of  tin- 
matrix.  Ultimately,  these  stress  concentrations  form  a  primary  source  of 
damage  initiation  and  property  degradation  particularly  in  the  presence  oi 
cyclic  loadings.  Recent  experimented  investigations  carried  out  by  Bal:i'  et. 
al  i  19SC)  on  graphite  epoxy  laminates  which  have  been  weakened  by  a  cir¬ 
cular  hole  give  us  a  better  insight  of  this  damage  growth  development  under 
the  action  of  cyclic  loadings.  In  general,  the  progression  of  this  damaged 
pn >cess  may  be  characterized  as  ( i )  debonding  along  fiber- matrix  interfaces 
(ii  i  matrix  cracking  parallel  to  the  fibers  (iii)  matrix  cracking  between  fibers 
(iv)  delamination  along  the  interface  of  two  adjacent  laminae  with  different 
fiber  orientations  and  (v)  fiber  breakage'. 

Thus,  if  rational  designs  using  hirer  reinforced -resin  matrix  composite 
laminates  are  to  be  made,  their  performance  under  static,  dynamic,  fatigue 
and  environmental  loads  need  to  bo  predieatable.  The  first  step  towards 
this  goal  is  the  realization  that  the  ultimate  failure,  as  well  as  many  other 
aspect s  of  the  composite  behavior,  is  the  result  of  the  growth  and  aeeumu- 
lation  of  micro-damage  to  the  fibers,  matrix  and  their  interfaces.  Thus, 
it  appears  that  any  generally  successful  model  of  performance  and  failure 
must  incorporate  the  effects  of  this  damage  in  some  way.  This  certainly 
represents  a  challenge. 

Delaminat  ion  has  long  been  recognized  as  one  of  the'  most  import  tint 
failure  modes  in  laminated  composite  structures.  The  growth  of  a  delam¬ 
ination  may  result  in  a  substantial  reduction  of  strength  and  stiffness  of 
the  laminate.  The  identification,  therefore,  of  such  locations  in  a  compos¬ 
ite  structure  is  of  great  interest  to  the  designer.  Experimental  studies  by 
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Pipes  et.  al.  (1973)  have  shown  that  the  delamiiiation  inode  of  failure  is 
most  likely  to  initiate  at  the  free  edges.  One  conjectures,  therefore,  that  t h< 
stresses  a t  the  intersection  between  a  free  edge  and  an  interface  may  well 
be  singular.  Indeed,  recent  analytical  investigations  (Wang  et.  ah.  1952. 
and  Zwiers  e.  al.  19S2)  on  straight  free  edges  show  that  a  stress  singularity 
exists  there  for  certain  types  of  laminates. 

Alternatively,  a  curved  free  edge  is  inherently  a  three-dimensional  prob¬ 
lem  which  presents  greater  mathematical  difficulties.  For  this  reason,  past 
analyses  have  been  based  primarily  on  finite  element  methods  with  standard 
finite  elements  (Raju  et.  al.  19S2)  as  well  as  elements  which  incorporate 
the  stress  singularity  in  the  formulation  (Rybicki  et.  al.  1975.  Erivson 
et.  al.  (1954).  While  such  methods  can  provide  us  with  stress  trends  in 
the  boundary  laryer  region,  it  is  rather  difficult  to  extract  from  them  with 
certainty  the  order  of  the  prevailing  stress  singularity  which  is  present  at 
the  material  interface.  Moreover,  experimental  investigations  carried  out 
on  straight  edges  (Pagano,  1974)  show  that  the  laminate  stacking  sequence 
can  effect  the  static  strength  of  the  laminates.  Similar  experimental  ob¬ 
servations  were  also  made  by  Daniel  et..  al  (1974)  on  plates  with  circular 
holt's.  The  subject,  therefore,  does  warrant  further  investigation. 

Recently,  the  author  investigated  analytically  the  interliminar  stresses 
at  the  boundary  layer  of  a  hole  free-edge  but  for  two  isotropic  materials 
of  different  material  constants  (Folias  19SS).  The  analysis  showed  that  the 
stress  field  there  possesses  a  weak1  singularity,  which  singularity  depends 
only  on  the  material  properties.  In  this  paper,  the  author  extends  this  anal¬ 
ysis  to  also  include  transversely  isotropic  laminae  with  a  [0'v / 90'  ] .  [0‘'/4o  j 
as  well  as  other  stacking  sequences. 

2  Formulation  of  the  Problem 

Consider  the  equilibrium  of  a  laminated  composite  plate  winch  occupies 
the  space  \.r\  <  cc,  |//|  <  oc  and  |c|  <  2 h  and  contains  a  cylindrical  hole  of 
radius  a  whose  generators  are  perpendicular  to  the  bounding  planes,  namely 
r  =  ±2 /;.  The  plate  consists  of  laminae  made  of  transversely  isotropic 
material  with  a  0° /90° /DO"' / 0°  stacking  sequence.  Let  the  plate  be  subjected 

‘loss  than  0.33 


i 
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to  a  uniform  tensile  load  a u  along  the  y-i ixis  and  parallel  to  the  hounding 
planes  (see  Fig.  1). 

In  the  absence  of  body  forces,  the  coupled  differential  equations  govern¬ 
ing  the  displacement  functions  u,  v.  and  w  are 


C  +c  —  +  (C  +c  \ 

11  d.r2  +Cw,(9y2  +Co5(9r2  +  (Cl2  +  C'*>)  Qj-dy 

+  (  3  +  C55  )  “  —  =  0 

0x0: 

r,  02u  ^  02r  ^  02  r  02r 

(C21  +  C  C.,;  )  0  +Cc,c,-r— 7  +  C22TT:  ~  C-hTT  7 

oxoij  ox-  oij-  a:- 

+  (^•23  +  £  1 1  )t7  7~  =  0 

OyOz 

r,  r  d’u  d2  r  d~ir 

(  C  31  +  Car, )  77  77“  +  (  C 32  +  C  )  —  -  r  C  57  “7  'T 

OXO  Z  OXOZ  OX- 


,,  d-W  ,,  0‘X 

+  £•11  -q—T  +  C 33^~7  —  0. 
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where  the  C',/s  are  the  material  constants  defining  a  layer  which  has  its 
libers  running  parallel  to  the  .r-axis.  For  the  next  layer,  the  fibers  will  be 
running  parallel  to  the  y- axis  and  the  governing  equations  will  be  obtained 
from  the  above  by  simply  interchanging  the  appropriate  coordinates. 

The  stress-displacement  relations  for  the  latyer  are  given  by  the  consti¬ 
tutive  relations 


As  to  boundary  conditions,  we  require  that: 


(4) 


at  r  -  ±2 h  :  the  surface  stresses  must  vanish 


(0) 


at  ;  =  i/(  :  the  displacements  and  suiface  stresses  must  match  <C' 

at  r  =  a  :  the  surface  stresses  must  vanish.  (  7 

Finally,  in  order  to  complete  the  formulation  of  the  problem,  the  loading 
conditions  far  away  from  the  hole  must  be  satisfied. 


3  Asymptotic  Solution  at  the  Interface 

I  lie  main  objective  of  this  analysis  is  to  derive  an  asymptotic  solution 
for  the  3D  stress  field  in  the  immediate  vicinity  of  the  region  where  the 
interface  between  two  laminae  meets  the  free  of- stress  surface  of  the  hole. 
Thus,  guided  by  a  general  analytical  solution  for  the  equilibrium  of  a  linear 
elastic  isotropic  layer  which  the  author  has  recently  constructed,  we  assume 
the  complementary  displacement  field  to  be  of  the  form 

( i )  for  lamina  [O'  ] 


//(1)  =  sin  0  <  C i 
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where 
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4 

1  15) 

( ;v2  =  C,,(Cn  sin*  8  4  Co;  cos*  W )  -f  CV,r.(  G„.  sin*  8  4  Cjj  cos*  8't  (  1C  ) 

t  ;t:j  =  C  1 1  C  .  lli' 

In  writing  the  above  displacement s.  we  used  a  cylindrical  coordinate 
system  and.  furthermore,  assumed  that  (r  —  a)  <<  a.  In  view  of 
the  .above,  one,  by  direct  substitution,  can  show  that  the  governing 
equations  ( 1 )— (3)  arc  indeed  satisfied  provided  the  unknown  function 
H  satisfies  the  differential  relation 
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whe:  c  the  fj.  f j  and  r:!  arc  functions  of  C\t  and  4  and  rc]>r<  lie 
root.-  of  the  cubic  equation 
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(  :i'a  —  C  -i-i  C  .  I  33  > 

and  the  function  H(r  —  a.h  —  : )  is  of  the  same  form  as  tlie  H  of  layer 
[O'] .  except  that  (q,  e->,  f:i  are  now  replaced  by  the  appropriate  q.  <y, 
e'i  of  layer  [90  ']. 

It  remains,  therefore,  for  us  to  construct  a  solution  to  equation  (IS).  To 
accomplish  this,  we  introduce  the  local,  to  the  corner,  stretched  coordinate 
system  (see  Fig.  2),  i.e. 


a  —  p  cc  >s  o 


(h  -  =) 


l>  Sill  o. 
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( hnitting  the  long  and  tedious  matliemataical  details,  the  solution  to  equa¬ 
tion  (IS),  in  terms  of  the  local  coordinate:-'.  is  found  to  be 
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where 
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in* i  n.  .4,  and  J3,(/  —  1.2.3)  arc  constants  to  he  detennined  from  the 
n  inndary  conditions. 

Sul  ist  it  ut  ing  t  he  previously  constructed  displacement  held  into  t  he  hound - 
: i  v  et  >ndit  ions: 


at  o  =  (): 


r(,»  =  t(2) 

tz  1  rz 

r(,)  =  r(2) 

y~  v: 


i  42 


43 


(44- 


( 40 ) 


(4G  i 
(47 ) 
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ct!'.)  =  o 


( 4: 


(i) 


i  40 


rll'  =0 


I  50 ' 


at  °=  o; 


=  o 


51 


_  ,, 

>r0  —  0 


=  0 


53 . 


we  arrive  ;it  a  system  of  twelve'  algebraic  equations.  the  determinant  of 
which  must  vanish.  This  latter  condition  leads  to  the  determination  of  the 
I'ha.racterist  ic  values  o.  In  General,  the  values  of  a  depend  on  the  material 


eonsfnats  C\j.  as  well  as  on  the  angle  H. 


4  The  Isotropic  Case 

A-  a  limit  check,  we  let  the  laminae  he  homogeneous  and  isotropic  hut  of 
dili'erent  material  constants.  Without  going  into  the  mathematical  details, 
the  material  constants  of  lamina  1  become: 


„  =  c,  =  c,  =  2(^;;)6', 

1 54  - 

C  l  =  c  55  —  Coe,  =  Ct 

(55 i 

2n, 

C 12  —  C 23  —  C\3  —  G\ , 

1  “  1V\ 

(5G) 

in  view  of  which  the  displacement  reduce  to: 
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H  =  //'  {.4]  cos( no)  +  B\  sin(ao)  +  -4;  cos(a  —  2)o 

( CO  i 

-1-  sin(o  —  2  )o  +  .4:3  cos(a  —  4)c>  +  Bi  sin(o  — 4)o  }  +  (/',’rl). 

Similarly,  the  H  collapses  to  the  same  expression  except  that  the  constant." 
.4,  and  Bl  tire  replaced  by  .4,  and  B ,  respectively.  The  nunierictd  results 
for  this  case  lead  to  the  same  results  as  those  recently  reported  by  Folia> 
( 1  OSS ).  Fig.  3,  for  example,  depicts  typical  values  of  a. 


5  Characterization  of  the  Free— Edge  Stress 
Singularity 

Returning  next  to  the  algebraic  system  (42)-(53).  if  one  considers  the  case 
of  a  graphite/epoxy  layer,  with  coefficients  C\}  (Knight  19S2) 

'  20.G22S  1.03S1  1.03S1  0.0000  0.0000  0.0000  ' 

1.03S1  2.2301  1.2301  0.0000  0.0000  0.0000 

_  1.03S1  1.2301  2.2301  0.0000  0.0000  0.0000 

ClJ  ~  0.0000  0.0000  0.0000  0.5000  0.0000  0.0000  '  } 

0.0000  0.0000  0.0000  0.0000  0.8690  0.0000 

0.0000  0.0000  0.0000  0.0000  0.0000  0.S69C  . 
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tln'ii  the  requirement  of  the  determinant  of  tin*  system  to  vanish  lead--  to  a 
1 ranscendent til  (‘([nation  for  the  roots  o.  The  only  roots  of  practical  interest 
are  those  which  lie  in  the  interval  5  <  ft  a  <  C.  The  numerical  results 
for  the  12  x  12  system  were  carried  out  in  double  precision.  Omitting 
the  long  and  tedious  numerical  details,  the  values  of  the  characteristic  o 
for  [0'790'j.  [0L'/70  ].  [0L'/45"]  and  [0’/20u]  interfaces*  are  shown  in  Fig. 
4.  Two  important  characteristics  are  worth  mentioning.  First,  the  stress 
singularity  is  a  funciton  of  the  material  constants  CtJ.  the  angle  of  sweep  9 
and  the  fiber  orientation  of  the  respective  laminae.  Second,  the  singularity 
strength  for  anisotropic  materials  appear  to  be  much  weaker  than  that 
of  isotropic  materials.  The  latter  nitty  have  severe  consequences  to  the 
damage  process  and  to  the  reduction  of  the  overall  strength  in  the  plate. 
As  a  practical  matter,  if  one  plots  the  max(u  —  G)  its  it  function  of  the  fibei 
untile  orientation  ,i  for  a  [0  ' / S:'\  interface  one  can  identify  the  most  and 
least  desirable  fiber  orientations.  This  is  depicted  in  Fig.  5. 

Similar  stress  singularity  profiles  (see  Fig.  4)  have  also  been  obtained  by 
\\  ang  et .  til.  ( 19S2)  in  their  pioneering  work  on  straight  edges  using  a  differ¬ 
ent  method  of  solution.  The  present  analysis  complements  that  and  shows 
that,  for  sufficiently  large  holes,  the  results  for  curved  edges  will  be  the 
"tune  as  those  obtained  near  it  straight  free  edge  provided  layer  orient  at  ion- 
in  the  second  problem  are  properly  chosen  to  reflect  the  circumpherential 
position  of  a  point  on  the  hole  boundary  and  the  interface.  While  this  re¬ 
mit  was  to  be  expected,  it  could  not  be  taken  for  granted.  This  is  because 
the  latter  method  represents  a  descrete.  rather  than  it  continuous,  approach 
and  the  outcome  of  the  limiting  process  laid  to  be  established.  Moreover, 
the  present  method  of  solution  shows  how  a  3D  analysis  can  indeed  be  de¬ 
veloped  to  also  include  this  continuous  dependence  on  the  angle  9  and  thus 
provide  us  with  further  insight  on  the  construction  of  such  3D  solutions  to 
transversely  isotropic  plates  with  more  complicated  Haw  geometries. 

It  may  further  be  noted  that  the  macromechanical  approach  actually 
underestimates  the  value  of  the  stresses  at  such  regions.  For  example,  if 
we  examine  the  local  geometry  from  a  micromechanical  point  of  view.  e.  g. 
at  9  =  (F  and  for  a  [0°/9UJ]  interface,  one  notices  that  the  adjacent  fiber 
of  layer  [90']  intersects  the  free  surface  of  the  hole  boundary  perpendicu¬ 
larly.  The  explicit  3D  solution  for  the  stress  field  in  such  regions  is  also 

‘"I  lic  results  are  presently  generalized  for  an  arbitrary  [0/d]  interface. 
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known.  In  part icular.  for  a  glass  fiber  embedded  into  an  epoxy  matrix  the 
stress  singularity  is  found  to  In'  0.24S9  ( Folias  19S9.  Fob  as  1990a)  whip 
for  a  carbon  fiber  embedded  into  an  epoxy  matrix  is  found  to  be  0.2S15 
(Folias  1990b).  The  former  analysis  assumes  the  fiber  to  be  of  an  isotropic 
material  while  the  latter  assumes  the  fiber  to  Ire  of  a  transversely  isotropic 
material.  Comparing  these  results  with  those  of  Fig.  4  it  is  clear  that  th< 
stress  singularity  predicted  by  the  macromechanical  theory  is  indeed  under¬ 
estimated.  Such  information  becomes  essential  for  the  proper  estimation  of 
the  local  damaged  zone.  This  matter  will  be  discussed  further  later  on. 

In  the  ease  of  a  [O' /90‘  ]  interface,  the  profile  of  the  characteristic  value 
o  versus  B  is  symmetric  with  respect  to  the  line  B  =  45''.  The  same  behavior 
was  also  obtained  by  Ericson  et.  al.  (19S4)  by  using  finite  elements.  T lie 
present  results,  however,  exhibit  a  stronger  singularity  than  that  found 
in  the  above  reference.  The  author  attributes  this  to  two  factors.  First 
the  material  constants  were  different  and  second  it  is  rather  difficult  to 
obtain  accurate  results  for  the  singularity  strength  based  on  finite  element 
analyses.  On  the  other  hand,  it  is  impressive  indeed  that  Ericson  et.  al. 
(19S4)  as  well  as  Raju  et.  al  \19S2)  were  able  to  recover  the  exact  profile 
as  a  function  of  B  and  the  relative  magnitude. 

In  order  to  make  a  proper  comparision  with  the  results  of  Ericson  et.  ah. 
one  should  use  the  same  material  constants.  CtJ .  as  those  used  in  the  above 
reference.  Computing,  therefore,  the  Cj/.s  from  the  data  of  the  above  ref¬ 
erence  (see  Appendix),  our  analysis  gives  the  characteristic  values  depicted 
in  Fig.  G.  At  B  =  45T  for  example,  a  =  5.9755  or  o  —  G  =  0.0245.  If  we 
now  compare  t h i .-»  value  with  that  found  by  Wanget.  al.  (19S2).  for  a  rfc  45 
straight  edge  interface,  i.e.  o  —  G  =  0.0255,  we  see  that,  the  comparison 
is  very  good.  The  minor  difference  is  probably  due  to  the  small  variation 
of  the  C\j  s  used  depending  as  to  how  they  are  computed.  Our  results  are 
based  on  the  C,/s  shown  in  the  Appendix.  The  results  in  the  region  be¬ 
tween  20  <  6  <  SO  compare  very  well  with  those  reported  by  Ericson  et. 
al.  On  the  other  hand,  for  0  <  6  <  20,  our  singularity  strength  is  found  to 
be  slightly  higher  and  t lie  characteristic  bell  shape  profile  is  preserved. 

Pagano  et.  al.  (1973)  have  shown  that  high  tensile  <7,-  stresses  are 
associated  with  the  decreased  laminate  strengths  reported  by  Pipes  et.  al. 
(1973).  This  observation  points  to  the  importance  of  understanding  the 
interlaminar  stress  behavior  near  free  edges  in  laminates.  It  is  now  possi¬ 
ble  to  compute  the  interlaminar  stresses  adjacent  to  the  hole  surface.  In 
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particular. 


rr.u  ~  !>'  f'M  (#.  o:  C\,-.  C,„„)  .  I  02  i 

where  the  /•>,  .are  rathem  long  and  complicated  functions  of  the  .angles  H 
and  o  and  the  material  constants.  Tin1  plots  for  a...  ~r: .  ri<:.  for  a  [()  '90 
inteiface,  are  given  in  Fig.  7.  Notice  that  the  max  a.,  stress  for  a  [0  /DO 
interface  occurs1  at  9  —  23'.  The  result  is  in  agreement  with  that  found 
previously  based  on  finite  (dements  (Ericson  et.  al  .and  Raju  et.  ah). 

The  re;td('r  may  also  note  that  .all  the  stresses  are  normalized  with  re¬ 
spect  to  a  function  C'(9)  which  is  negative  and  may  vary  slightly  as  one 
travels  along  the  direction  of  9.  This  function  of  theta  is  to  be  determined 
from  a  separate  analysis  valid  across  the  thickness  of  the  plate  by  comparing 
the  relative  a,-,,,  stress  functions.  Tin*  analysis  is  presently  being  completed 
and  will  be  reported  separately.  The  present  analysis,  however,  does  pro¬ 
vide  us  with  the  relative'  magnitude's  of  the  interlaminar  stresses  t  see  Fin. 

Si- 

The  results  show  the  shear  stress  rT:  to  be  relatively  low  throughout. 
On  the  other  hand,  the  normal  sfess  a::  appeals  to  be  dominant  in  the 
range  —30"  <  6  <  30"  beyond  which  the  shear  stress  t,-,~  now  becomes  the 
controlling  stress,  i.e.  along  30"  <  9  <  GO"  with  its  maximum  occuring  tit 
9  —  45'  .  A  more  clear  plot  for  the  stress  a::  at  the  interface  is  given  in 
Fig.  9.  Finally.  Fig.  10  depicts  the  jump  which  exists  on  the  stress  er.-..,  as 
one  moves  across  the  interface,  i.e.  at  o  =  0+  and  O  —  0_.  The  maximum 
difference  occurs  at  9  —  0  and  is  approximately  27 c/e.  In  Figs.  11  and  12 
we  compare  the  stress  erf-  at  6  =  0—  and  6  =-90"  and  at  o  =  0+  and 
O  =  +90  respectively. 

In  view  of  the  above,  one  may  draw  the  following  conclusions  for  a 
[0/90"]  stacking  sequence: 

(a)  In  the  vicinity  of  the  interface  there  exists  a  boundary  laryer  effect 
where  the  stress  field  changes  rather  abruptly. 

(b)  The  res  k  of  delamination  initiation  is  highest  at  9  equal  to  23". 

'*  I’lie'  rcaele'r  should  also  take’  into  account  the  fact  that  (6- a)  is  a  maximum  at  0  = 
and  that  C(0)  is  negative.  The  results  are  based  on  data  give  b\  Eq.  (61).  Also  at  0  —  0, 
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(c)  As  one  mows  away  from  the  plane  of  the  interface,  the  stress  concen¬ 
tration  factor  in  layer  [90  ]  decreases  rapidly  (see  region  —20  <  ^  *  20 
where  debonding  along  fiber  matrix  interfaces  are  most  likely  to  ini¬ 
tiate). 

( <1 )  Substantial  damage  is  expected  to  occur  in  the  region —40  '  <  ft  <  40  . 

(e)  For  the  given  set  of  material  const  tints  C',,.  n  -  G  is  maximum  .at 
0*  =0;,.  =td/2). 

(f)  In  general,  the  magnitude  of  the  signularity  strength  depends  on  the 
material  constants.  Cu.  and  on  t lie  fiber  orientation  of  the  two  re¬ 
spect  ive  laminae. 

(g)  As  one  moves  ;i]>]>roximately  one  radius  distance  awtiy  from  the  hole 
surftice,  the  s.c.f.  (stress  concent  art  ion  ftictor)  is  expected  to  decrease 
to  within  lO'/f  of  the  value  of  the  corresponding  case  of  a  plate  wit  bout 
a  hole  (Folias  1990  isotropic  case). 


6  Discussion 

Delamination  at  free- edges  has  been  a  problem  of  significant  concern  in 
the  design  of  advanced  fiber  composites.  The  separation  of  the  laminae, 
caused  by  high  local  interlaminar  stresses  and  low  strength  along  the  lam¬ 
inae  interface,  nitty  result  in  ineffective  load  transfer,  reduction  in  stiffness 
and  ultimately  loss  of  structural  integrity.  In  this  study,  the  problem  has 
been  investigated  by  treating  each  lamina  as  a  homogeneous,  transversely 
isotropic,  material.  Thus  the  micromechanics  effects  of  fiber  size  are  not 
included,  although  a  few  of  these  effects,  eg.  when  a  fiber  meets  a  free 
surftice  have  been  reported  separately  (Folias,  1989). 

The  analytical  investigation  of  the  3D  stress  field  adjacent  to  the  hole 
and  in  the  vicinity  of  the  interfaces  of  two  laminae  of  [0u/90°],  and  [0U/45U], 
and  other  fiber  orientation  shows  t lie  stresses  to  be  singular.  atJ  ~  //Wl.  In 
general,  the  singularity  exponent  depends  on  the  material  properties,  the 
corresponding  fiber  orientations  as  well  as  the  angle  of  sweep.  The  results 
provide  us  with  better  insight  for  the  proper  understanding  of  interlami¬ 
nar  stresses  and  the  effect  which  they  have  on  the  mechanism  of  failure. 
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For  example,  for  isotropic  laminae  the  stress  singularity  is  a  rather  weak 
singularity  (<  0.33),  while  in  transversely  isotropic  laminae  it  can  he  of  a 
much  higher  order  (<  0.99)  suggesting,  therefore,  a  greater  influence  in  the 
damage  process. 

Pagano  et.  al  (1971,  1974)  postulated  that  throughout  the  thickness  it 
is  the  stress  rr::  which  is  the  main  cause  of  delamination  for  polymeric  based 
structural  composites.  This  observation  w;is  based  on  experimental  its  well 
its  analytical  evidence  for  laminates  with  different  stacking  sequences.  Our 
analytical  results  are  consistent  with  this  observtion  particularly  in  the 
region  where  0  is  small.  Alternatively,  for  larger  values  of  0  the  shear  stress 
Tfi~  plays  a  much  more  dominant  role  on  interlaminar  failures.  Moreover, 
in  the  case  of  a  [0  /90  ]  interface  delamination  is  most  likely  to  take  place 
at  0  —  23'.  The  strain  energy  release  rate  may  now  be  used  in  conjunction 
with  the  local  stress  field  to  predict  delamination  failures.  This  matter  is 
presently  under  investigation. 

Finally,  one  may  conclude  that  it  is  possible  to  reduce  the  likelyhood 
of  the  delamination  mode  of  failure  and  thereby  increase.1  the  laminate 
strength.  This  can  be  accomplished  by  carefully  considering  the  effects  of 
the  singularity  curves,  the  stress  curves,  the  load  direction  and  the  individ¬ 
ual  fiber  orientations  at  each  interface'.  Mon-over,  in  future  applications  it 
may  be  possible  to  choose  the  material  constants  C'tJ  so  that  tlu-  coefficients 
of  the  singular  terms  of  the  interlaminar  stresses  vanish. 

7  Micromechanical  Considerations 

If  one  examiners  the  region  adjacent  to  the  hole  and  at  the  vicinity  of  the1 
interface  from  a  micromechanical  point  of  view  one  notices  the  arrangement 
depicted  in  Fig.  13(a).  More  specifically,  in  one  layer  the  fibers  intersect 
the  free  surface  of  the  hole  perpendicularly,  while  in  the  other  layer  they  run 
tangent  to  its  surface  (see  Fig.  13b).  Locally,  however,  this  represents  the 
problem  of  a  fiber  intersecting  a  free  surface  and  being  subjected  to  a  biaxial 
load  normal  to  the  direction  of  the  fiber.  As  part  of  this  investigation,  this 
problem  has  been  recently  studied  by  the  author  and  the  results  for  the  case 
of  an  isotropic  fiber  have  been  reported  in  Folias  (19S9).  The  explicit  local 
stress  field  for  one  fiber  may  be  found  in  part  B  of  this  report.  Without 
going  nto  the  mathematical  details,  the  local  stress  field  was  found  to  be 
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singular.  More  specifically,  the  stresses  were  shown  to  be  of  the  form 


(7xj  ~  f>("  (03  i 

where  the  exponent  d  in  general  is  a  function  of  the  material  properties.  In 
the  case  of  a  glass  fiber  embedded  into  an  epoxy  matrix  with  the  properties 

Gj  =  35.006'Pu  v;  =  0.22 
G,„  =  2.10GPa  um  =0.34 

the  stress  singularity  strength  is  found  to  be  equal  to  2  —  o  =  0.24S9.  It  is 
well  recognized,  however,  that  carbon  fibers  are  transversely  isotropic  rather 
than  isotropic.  For  this  reason,  the  previous  analysis  has  been  extended  to 
also  include  a  transversely  isotropic  fiber  (Folias  and  Li  1990).  Without 
going  into  the  mathematical  details,  the  stress  singularity  for  a  carbon 
fiber  emb  edded  into  an  epoxy  matrix  was  found  to  be  slightly  higher,  i.e. 

2  -  d  =  0.2S2. 

Comparing  these  results  with  those  of  Fig.  4  (see  Fig.  14)  it  becomes 
clear  that  the  derived  stress  field  adjacent  to  the  hole  based  on  niacronie- 
chanical  considerations  is  indeed  underestimated.  Such  information  is  es¬ 
sential  if  one  is  to  predict  failure  initiation  and  local  damage  degradation 
adjacent,  to  the  hole. 

Thus  one  may  view  the  process  of  delamination  as  the  interfacial  crack¬ 
ing  of  the  last  row  of  fibers  which  intersect  the  free  surface  of  the  hole  (see 
Fig.  13b  which  corresponds  to  8  —  90', or  8  =  0  i£  the  figure  in  flipped). 
Recalling  the  results  of  the  surface  displacement  of  the  problem  studied  by 
Penado  and  Folias  (19S9,  sec  Figs.  14  and  15  of  that  paper),  one  concludes 
that  debonding  is  most  likely  to  take  place  at  the  lower  side  of  the  fiber, 
e.g.  see  Fig.  13b.  It  would  be  of  great  practical  interest  if  one  could  pre¬ 
dict  a  priori  the  critical  applied  stress  a0  which  may  cause  a  fiber/niatrix 
debonding  to  initiate  and  subsequently  propagate  and  become  a  full  size 
delamination  crack  (see  Fig.  7  of  Part  II).  In  order  to  accomplish  this,  we 
use  the  results  of  part  II,  in  particular  equations  (29b  and  33)  and  apply 
them  at  the  position'  8  %  0.  Thus 

‘’Notice  that  a.z  is  ss  the  same  for  0  <  0  <  25. 
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( 0 .0944  1(9.595$  )(44rru  ){  ~  r^Joaj.i 

It  V 

(  C-t  ! 

—  l.Sl-fGyrfu’G,,,  • 

where  the  first  term  is  due  to  the  periodicity  of  the  fibers,  the  third  term  is 
due  to  the  local  stress  concentration  factor,  a  is  the  radius  of  the  hole.*//  is 
the  thickness  of  the  layer,  (ij  the  radius  of  the  fiber.  ->12  is  the1  surface  enemy, 
and  G,„  the  shear  modulus  of  the  matrix  and  2 .i  the  angle  of  debond.  Xexr. 
solving  for  the  critical  stress  au  one  finds 

(<7())rr  ~  C.1S9  X  10  t  /  —  —  [  ~  )  .  (Go  1 

V  2 aj.i  a 

Assuming  a  fiber  diameter  of  0.001cm.  a  matrix  shear  modulus  of  2.1  GPu. 
and  a  surface  energy  of  270 ,J /nr  for  well-bonded  interfaces,  one  finds 

(<7ubr  -  2.00— t==(  —  f)A' M Fa.  (CGI 

v/2  i  a 

For  Poorly  bonded  interfaces  the  constant  is  approximately  one  half  of 
that. 
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Appendix 


From  Eriesoii  et.  al  data  (1984) 

En  =  13SGTF;  =  E:„  =  14.567’./ 

Gj  j  =  6 13  —  6  =  5.9  Ci  Pa 


i'u  ~  7'ih  —  t'j  f  =  9.91 . 


we  compute: 


6]  i  =  139.0351 

C,  2  =  3.9002 

Ci.j  =  3.9002 

Cjj  =  15.2779 

Cj .  =  3.2944 

Get  =  lo.2i  <9 

C. I  I  =  C  =  G.r.  =  0.9 
in  view  of  which  our  program  then  gives: 


0 

a 

0.1 

5.9440 

10 

5. 9409 

15 

5.909S 

20 

5.9714 

30 

5.9739 

45 

5.9755 
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FIGURE  CAPTIONS 


Fig.  1.  Laminated  plate  of  arbitrary  thickness  with  a  circular  hole. 

Fig.  2.  Definition  of  local  coordinates  at  the  interface. 

Fig.  3.  Singularity  strength  for  isotropic  laminae. 

Fig.  4.  Singularity  strength  for  transversely  isotropic  laminae  [0 '  /90  ] . 

Fig.  5.  Best  and  worse  fiber  orientations. 

Fig.  G.  Singularity  strength  for  transversely  isotropic  laminae  T)  /90  baset 
on  Ericson  et.  al.  data. 

Fig.  7.  Interlaminar  stresses  for  a  [0  /90  j  interface. 

Fig.  S.  Interlaminar  shear  stress  ratio  for  a  [0”/90]  interface. 

Fig.  9.  The  interlaminar  stress  a:z  vs  0. 

Fig.  10.  The  stress  concentration  factor  for  o  =  0+  and  0"  . 

Fig.  11.  The  stress  crjlj  for  6  —  0+  and  o  -  90”. 

Fig.  12.  The  stress  cr\'^  for  6  —  0“  and  o  =  -90  . 

Fig.  i3.  (a)  Geometrical  representation. 

(b)  Local  geometry  based  on  micromechanical  considerations.  { Po 
sit  ion  0  =  90”) 

Fig.  14.  Comparison  of  stress  singularities. 


20 


z 


°° 

90° 

90° 

0° 

Fig.  1.  Laminated  plate  of  arbitrary  thickness  with  a  circular  hole. 


FRONT  SURFACE 


Fig.  2.  Definition 


of  local  coordinates  at  the  interface 


Fig.  13.  (a)  Geometrical  representation 


MICROMECHANICAL  CONSIDERATIONS 


hole  surface 


Fjrr  ;3  (b)  Local  geometry  based  on  micromechanical  considerations.  (Po¬ 
sition  6  —  90°) 


Fig.  14.  Comparison  of  stress  singularities. 


Part  II 


ON  THE  PREDICTION  OF  FAILURE  AT  A  FIBER/ MATRIX  INTERFACE  IN 
A  COMPOSITE  SUBJECTED  TO  A  TRANSVERSE  TENSILE  LOAD 


ABSTRACT 


This  paper  deals  with  the  3D  stress  field  of  a  cylindrical  fiber  which  is 
embedded  into  a  resin  matrix.  The  composite  is  then  subjected  to  a  uniform 
tensile  load  ao-  The  strain  energy  release  rate  is  computed  and  the  criterion  is 
used  to  predict  debonding  initiation  at  the  fiber/matrix  interface.  The 
analysis  shows  that  this  failure  is  most  likely  to  occur  at  the  free  surface,  ie  the 
region  where  the  fiber  intersects  a  free  surface  for  example  a  hole,  an  edge,  or 
a  crack.  Moreover,  it  will  occur  at  approximately  (1/10)  the  load  value 
required  for  the  same  failure  to  commense  at  the  center  of  the  fiber  length. 

The  results  are  also  extended  to  include  a  doubly  periodic  array  of  fibers 
which  are  embedded  into  a  matrix.  Based  on  3D  considerations,  the  stiffness 
matrix  is  shown  to  increase  as  the  volume  fraction  of  the  fibers  increases. 
Similarly,  the  stress  arr  in  the  matrix  is  shown  to  decrease  as  the  volume 
fraction  of  the  fibers  increases. 
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INTRODUCTION1 


It  is  well  recognized  that  fiber  composite  materials  are  verv  attractive 
for  use  in  aerospace,  automotive  and  other  applications.  These  composites 
consist  of  relatively  stiff  fibers  which  are  embedded  into  a  lower  stiffness 
matrix.  Although  in  most  designs  the  fibers  are  aligned  so  that  they  are 
parallel  to  the  direction  of  the  external  loads,  it  is  almost  impossible  to  avoid 
induced  transverse  stresses  which  may  lead  to  premature  failure  of  the 
laminate.  An  excellent  example  of  this  is  the  case  of  a  filament  wound 
pressure  vessel  in  which  the  presence  of  a  curvature  induces  bending  as  well 
as  transverse  stresses  (Folias,  1965).  However,  in  order  to  be  able  to  predict 
their  failing  characteristics,  particularly  in  the  neighborhood  of  free  surfaces 
such  as  holes,  edges  etc.,  it  is  necessary  to  know  the  local  stress  behavior  from 
a  3D  point  of  view. 

An  overall  summary  of  some  of  the  results,  which  are  based  on  2D 
elasticity  considerations  can  be  found  in  the  books  by  Hull  (1981)  and  by 
Chamis  (1975).  In  their  pioneering  work,  Adams  and  Doner  (1967)  used  finite 
differences  to  solve  the  problem  of  a  doubly  periodic  array  of  elastic  fibers 
contained  in  an  elastic  matrix  and  subjected  to  a  transverse  load.  Their 
results  reveal  the  dependence  of  the  maximum  principal  stress  versus  the 
constituent  stiffness  ratio  (Ef/Em)  for  various  fiber  volume  ratios.  A  few 
years  later,  Yu  and  Sendeckyj  (1974)  used  a  complex  variable  approach  to 
solve  the  problem  of  multiple  inclusions  embedded  into  an  infinite  matrix. 
Their  results  were  subsequently  specialized  to  cases  of  two  and  three 
inclusions  thus  providing  us  with  futher  insight  into  the  strength  of  the 
composite.  On  the  other  hand,  the  separation  of  a  smooth  circular  inclusion 
from  a  matrix  was  investigated  by  Keer,  Dundurs  and  Kiattikomol  (1973).  By 
using  finite  integral  transforms,  they  were  able  to  reduce  the  problem  to  that 
of  a  Fredholm  integral  equation  with  a  weakly  singular  kernel.  Thus, 
extracting  the  singular  part  of  the  solution,  they  were  able  to  reduce  the 
remaining  problem  to  a  simpler  one  which  lends  itself  to  an  effective 
numerical  solution.  Their  results  are  very  general  and  are  applicable  to 
various  combinations  of  material  properties  and  loads. 
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In  this  paper,  use  of  the  local,  3D,  stress  field  will  be  made  in  order  to 
examine  the  dependence  of  the  stress  crr,  in  the  matrix,  on  the  ratio  (Gf/Gm). 
The  strain  energy  release  rate  will  then  be  computed  in  order  to  predict  crack 
initiation  at  the  fiber/matrix  interface.  Particular  emphasis  will  be  placed  in 
the  region  where  fibers  meet  a  free  surface  as  well  as  at  the  center  of  a  fiber's 
length. 

FORMULATION  OF  THE  PROBLEM 


Let  us  consider  a  cylindrical  fiber  of  homogeneous  and  isotropic 
material,  e.g.  a  glass  fiber,  which  is  embedded  into  a  matrix  of  also 
homogeneous  and  isotropic  material. 

Futhermore,  we  assume  the  matrix  to  be  a  rectangular  plate  with  finite 

dimensions  2w,  2/,  and  2h  as  defined  by  fig.  1.  For  simplicity,  we  assume 
w  / 

—  >  8  and  -  >  8.  Such  an  assumption  will  guarantee  that  the  boundary  planes 

x  =  ±  w,  and  y  =  ±  /,  will  not  effect  the  local  stress  field  adjacent  to  the  fiber.* 
Thus,  mathematically,  one  may  consider  the  boundaries  in  the  x  and  y 
directions  to  extend  to  infinity.  As  to  loading,  the  plate  is  subjected  to  a 
uniform  tensile  load  a0  in  the  direction  of  the  y-axis  and  parallel  to  the 
bounding  planes  (see  Fig.  1). 

In  the  absence  of  body  forces,  the  coupled  differential  equations 
governing  the  displacement  functions  up  are 


1 


1-2  Vi 


de<j) 

dx[ 


+  V2 


=  0,  i  =  1,2,3,  j  =  1,2, 


(1) 


where  V2  is  the  Laplacian  operator,  vj  is  Poisson's  ratio,  u|^  and  uj2^  represent 
the  displacement  functions  in  media  1  (matrix)  and  2  (fiber)  respectively,  and 


e(j) 


9  u^ 

"9xi~  ;  1  =  1,2,3 ;  j  =  1,2‘ 


(2) 


The  stress-displacement  relations  are  given  by  Hooke's  law  as 


*  This  can  be  seen  from  the  results  which  were  recently  reported  bv  (Pcnado  and  Folias  (1989). 
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(3) 


$  * x)  4k  5i/ + 2g(  •  4’- 

where  Aj  and  Gj  are  the  Lame'  constants  describing  media  1  and  2. 

THE  SOLUTION  FOR  ONE  FIBER 

A.  Region  where  fiber  intersects  the  free  edge 

This  problem  was  recently  investigated  by  the  author  (Folias  1989)  who 
was  able  to  recover,  explicitly,  the  three  dimensional  stress  field  adjacent  to 
the  surface  of  the  fiber’  .  Without  going  into  the  mathematical  details,  the 
displacement  and  stress  fields  for  the  matrix  are  given  in  terms  of  the  local 
coordinate  system  (see  fig.  2)  by: 

(i)  displacement  field: 

u(1)  =  An  pa*l  sinG  {b  [2  (1-vi)  cos  (oc-l)(J>  -  (a-1)  sino  sin(a-2)o]  (4) 

-  (a+1)  [(l-2vi)  sin  (a-l)0  +  (a-1)  sin0  cos  (a-2)  0]  ] cos(2nG) 

v(1)  =  An  Pa~^  cosG  {b  [2  (1-vj)  cos  (a- 1)0  -  (a-1)  sino  sin(a-2)0  ]  (5) 

-  (a+1)  [(l-2vi)  sin  (a-1)  0  +  (a-1)  sin0  cos  (a-2)0]  )  cos(2n0) 

w^)  =  An  pa'^  {b  [-  (l-2vi)  sin  (a-l)0  +  (a-1)  sin0  cos  (a-2)  0  ]  (6) 

-  (a+1)  [  2  (1-vi)  cos  (a-l)0  +  (a-1)  sin0  sin(a+2)  0  ]}cos(2n0) 

(ii)  stress  field: 

°rr  =  (a-D  Pa’^  {b  [  2  cos  (a-2)0  -  (a-2)  sinO  sin(a-3)p]  (7) 

-  (a+1)  [sin  (a-2)0  +  (a-2)  sin0  cos  (a-3)0  ]  }  cos(2n0) 

000  =  4v]G^^  (a-1)  An  pa'^  {b  cos  (a-2)0  -  (a+1)  sin(a-2)d)  cos(2n0)  (8) 


A  similar  analysis  for  a  transversely  isotropic  fiber  meeting  a  free  surface  has  recently  been 
completed  and  the  results  will  be  reported  soon. 


(9) 


<^2z  =  (a-1)  An  pa'2  {b  (a-2)  sinp  sin(a-3)0 

+  (a+1)  [(a-2)  sin0  cos  (a-3)0  -  sin(a-2)p  ]}  cos(2n9) 

x^  =  2G^  (a-1)  An  pa*^  {b  [sin  (a-2)<J>  +  (a-2)  sino  cos  (a-3)o]  (10) 
-  (a+1)  (a-2)  sin4>  sin(a-3)p  }  cos(2n0) 

(nwi2) 


where  n  =  0,1,2,...  and  B  is  a  function  of  the  material  constants  and  An  is  a 
constant  to  be  determined  from  the  boundary  conditions  far  away  from  the 
fiber*.  In  general,  the  characteristic  value  of  a  depends  on  the  material 
constants  of  the  fiber  as  well  as  of  the  matrix.  A  typical  example  is  given  in 
fig.  3. 

Upon  examination  of  the  stress  field,  the  following  remarks  are  worthy 
of  note.  First,  the  stress  field  in  the  neighborhood  where  the  fiber  meets  the 
free  surface  is  signular.  Moreover,  in  the  limiting  case  of  a  perfectly  rigid 
inclusion  this  singularity  strength  reaches  the  value  of  0.2888.  Second, 
boundary  conditions  and  are  satisfied  as  a  consequence  of  the 

odd  functional  behavior  in  0,  which  points  to  the  presence  of  a  boundary 
layer  solution  as  one  approaches  the  free  surface.  Third,  on  the  free  surface 
the  radial  stress  is  (1/vi)  times  the  circumpherential  stress.  This  suggests, 
therefore,  that  if  a  crack  was  to  initiate,  it  would  propagate  along,  (or  very 
very  close  to)  the  fiber/matrix  interface.  Clearly,  the  occurance  of  either 
adhesive  or  cohesive  failure  will  depend  on  the  relative  strengths  of  the 
interface,  of  the  fiber,  and  of  the  matrix.  All  things  being  equal,  the  analysis 
shows  the  stresses  to  be  highest  at  the  interface,  thus  pointing  to  an  adhesive 
type  of  failure. 


*  for  one  fiber  n=0,l.,  while  for  a  periodic  extension  n=0,l,2,.. 
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B.  Interior  region 

The,  3D,  stress  field  for  this  region  has  also  been  recovered  by  Penado 

and  Folias  (1989)  and  the  results  for  various  (a/h)  and  (G?/G])  ratios  may  be 

found  in  the  literature.  The  results  have  subsequently  been  extended  (Folias 

and  Liu,  1990)  to  also  include  a  layer  of  modified  matrix  around  the  fiber. 

Thus  for  vi  =  0.34,  V2  =  0.22  and  (G2/G1)  =  16.67  the  stresses  o^  and  oil?  at  r  = 

1  r  uu 

a  and  for  all  Izl  <  h  are  given  in  figs  4  and  5  respectively.  Finally,  fig.  6  (for 

(1) 

X=0)  shows  the  variation  of  the  stress  crr  as  a  function  of  the  ratio  (Go/Gi). 

INTERFACE  FAILURE  CLOSE  TO  THE  FREE  SURFACE 

A  closer  inspection  of  the  local  stress  field  shows  that  a  crack  is  most 
likely  to  initiate  at  the  location  0  =  0  and  subsequently  propagate  along  the 
fiber/matrix  interface  until  it  reaches  a  nominal  value  of  the  arc  length 
beyond  which  it  will  advance  into  the  matrix.  Moreover,  once  the  crack 
begins  to  propagate,  it  will  simultaneously  propagate  along  the  interface  and 
parallel  to  the  axis  of  the  fiber  (mode  III).  Thus,  crack  propagation  will  be 
governed  initially  by  a  mode  I  failure  and  subsequently  by  a  combination  of 
mode  I  and  mode  III  failure.  It  is  now  possible  for  us  to  examine  the  first 
stage  of  the  failing  process  and  to  obtain  an  estimate  of  the  debonded  arc 
length  as  well  as  an  estimate  of  the  critical  transverse  stress  for  crack 
initiation. 

As  a  practical  matter,  we  will  consider  the  special  case  of  a  glass  fiber 
embedded  into  an  epoxy  matrix  with  the  following  properties 

Gi  =  2.10  GPa  vi  =0.34  (13) 

G2  =  35.00  GPa  v2  =  0.22  . 

Without  going  into  the  numerical  details,  the  constants  a,  A  and  B  for  this 
example  are  found  to  be* : 

*  The  constant  A  has  been  determined  by  comparing  the  displacement  w(1  *,  as  well  as  the 
(1) 

stress  crr  at  8  =  0,  at  z  =  h  and  for  (a/h)  =  0.5  with  the  work  of  Penado  and  Folias  (1989). 

A  s  X  An. 
n 
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a  =  1.7511,  G(D  A  a«'2  =  0.6349  a0  ,  B  =  2.1302  , 


(14) 


where  o0  now  has  the  units  of  GPa.  Thus,  from  equations  (7)-(12^  one  has 
(i)  at  9  =  0  and  9  =  0: 


41?  *  4.0633  c0 

^r'aV0.2489 

a 

(15) 

X  / 

(1)  (1) 
°99  =  V1  °rr  ' 

(16) 

tc/2  and  9  =  0: 

=  1.9163  o0 

f'h-z>0.2489 

lh  J 

(17) 

°ee '  °-4844  <£’ 

(18) 

-  0.2930  o^- 

(19) 

It  is  clear  now  from  equations  (15)  and  (16)  that  crack  failure  is  most 

likely  to  initiate  and  subsequently  propagate  along  the  fiber/matrix  interface 

rather  than  perpendicular  to  it.  Similarly,  equations  (17)  and  (19)  suggest  that 

failure  in  the  direction  parallel  to  the  axis  of  the  fiber  is  dominated  first  by  a 

mode  I  and  second  by  a  mode  III  type  of  failure.  It  may  also  be  noted  that 

(1) 

arr  attains  a  maximum  at  9  =  0  and  descreases  as  one  travels  along  the 
surface  of  the  fiber. 

Finally,  based  on  3D  considerations,  the  stress  field  away  from  the 
edges,  z  =  ±h,  and  in  the  interior  of  the  plate  was  shown  to  be  non-singular 
(Penado  and  Folias  1989,  Folias  and  Liu  1990)  with* 

=  0.4090  =  0.4090  (1.4281  o0)  =  0.5841  c0  ,  at  0  =  0  (20a) 


*  These  results  are  valid  for  a  ratio  of  (a/h)  =  0.05  and  subject  to  the  assumption  that  (—  )  >  8 

and  (//a)  >  8  in  which  case  the  end  boundaries  in  x  and  y  have  insignificant  effects  on  the 
local  to  the  fiber  stress  field. 
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at  r  =  a  and  z  =  0.  Comparing  this  value  with  that  of  the  corresponding  plane 
strain  solution 


Gqq  =  j—  a;- =0.5152  Go,  at0  =  O  (20b) 

one  notices  that  it  is  approximately  13%  higher  in  value  due  to  the  presence 
of  the  stresses  in  the  third  dimension. 

It  is  now  possible  for  us  to  obtain  an  approximate  criterion  for  de¬ 
bonding  along  the  fiber/matrix  interface  based  on  Griffith's  theory  of  fracture. 
Thus,  following  the  work  of  Toya  (1974),  if  one  assumes  the  presence  of  an 
interface  crack  of  length  2a(3  and  if  futhermore  takes  into  account  the  local  3D 
stress  field,  then  Toya's  result  may  be  written  as 

(1/16)  (1.1337o0  )2  k  a  Ai  (1+4  e2)  n  N0  N0sin  p  exp  [2  e  (jt-p)]  =  2y  (21) 
where 


l+k2 

k  ~  l+k2  +  (1+ki)  (G2/Gi) 


3-4vi 

hi  =  3~  vi 

1  +  vi 


for  plane  strain 
for  plane  stress 


Al 


k 

4Gi 


l+ki  +  (l+k2)  (G1/G2) 


(22) 

(23) 

(24) 


J_  /  [l+h2  (Gi/G2)  ] 
2k  L  ki  +(Gi/G2)  J 


(25) 


Nq  =  Go  '  k 


2  (1-k) 

k 


l+k2(G]/G2)  f  1 

kT^(Gi/G2)  exP  LP (2e  ‘  l}J 


Go  = 


1-  (cos  p  +  2e  sin  P)  exp  [2e  (n-p)]  +  (1-k)  (1  +  4e2)  sin2p 
2-k-k  (cos  P  +  2e  sinP)  exp  [2e  (71  -  P)] 


(26) 

(27) 
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where  No  is  the  complex  conjugate  of  N0/  y]2  is  the  specific  surface  energy  of 

the  interface  and  (3  the  angle  of  debonded  interface  (see  r!c,  while  it  is 
true  that  this  type  of  approach  does  not  provide  results  for  the  exact  initiation 
of  an  interface  crack  problem,  ie  from  a  condition  of  perfectly  bonded  interface 
to  that  of  a  partially  debonded  interface,  it  does,  however  provide  a  very  good 
first  approximation  to  this  complex  phenomenon.  The  author  is  well  aware 
of  that  and  is  presently  continuing  his  work  along  such  lines  and  with  some 
promise. 

Upon  rearranging,  equation  (21)  can  be  written  in  the  form* 

2y12 

~2~  =  (1.2853)  F  (Vi ,  Gi ;  (3),  (28) 

o0a 

where  F  is  a  function  of  the  material  constants  and  the  angle  |3  of  the 
debonded  interface.  A  plot  of  this  equation  for  conditions  of  plane  stress,  as 
well  as  of  plane  strain,  is  given  in  Fig.  8.  In  both  cases  the  maximum  occurs  at 
[3  =  60°.  Beyond  this  angle,  the  crack  will  gradually  curve  away  from  the 
interface  and  into  the  matrix. 

In  order  for  us  to  obtain  an  estimate  for  the  critical  stress  for  crack 
initiation  we  let  [3  0+,  ie  very  small  but  not  zero.  Thus,  for  our  example 

(a0)cr  V2  a  (3  =  1.8186  ^y^Gu;  at  z  =  0.  (29a) 

On  the  other  hand,  in  the  neighborhood  of  the  free  surface,  the  applied  stress 
is  much  higher  because  of  the  singularity  presence.  In  order  to  overcome  this 
difficulty,  one  may  average  the  local  stress  over  a  distance  equal  to  10%  of  that 
of  the  radius,  ie. 


*  It  should  be  noted  that  at  the  crack  ends  the  stress  field  oscillates  and  that  some  overlap  of 
the  crack  faces  takes  place.  This  matter  is  well  recognized  and  has  been  documented  by 
Williams  (1952),  Rice  et  al.  (1965)  and  England  (1965).  The  region  where  this  occurs,  however, 
is  so  small  (less  than  a  x  10‘3)  that  eq.  (28)  provides  a  good  approximation. 
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(°o)eff  "  (0.1  a) 


(30) 


°'la(4.0633  a0) 
JO 


\  -0.2489 


dq 


=  9.5958  c0  . 


Thus* 


(9.5958  o0)cr  V2  a  (3  =  1.8146  G]  ;  at  z  =  ±  h.  (29b) 

Combining  next  eqs.  (29a)  and  (29b)  one  finds 
(oo)cr  |at  z=h 

-  -  0.10,  (31) 

(°o)cr  | at  z=0 

ie.  the  critical  loading  stress  which  may  cause  failure  close  to  a  free  surface  is 
approximately  (1/10)  of  the  critical  stress  required  to  cause  the  same  failure  at 
the  center  of  the  fiber's  length.  Thus,  all  things  being  equal,  a  crack  will 
initiate  at  the  free  surface  and  will  propagate  along  the  periphery  of  the 
fiber/matrix  interface  as  well  as  parallel  to  the  axis  of  the  fiber. 

Focusing  next  our  attention  on  the  advancement  of  the  crack  along  the 
periphery  of  the  fiber  we  conclude  that  the  crack  will  advance  itself  to  a 
cricital  angle  of  p  =  60°.  Once  the  crack  has  reached*  9  =  60°,  the  local 
geometry  is  similar  to  that  of  a  hole.  This  problem  has  also  been  investigated 
for  the,  3D,  stress  field  close  to  a  free  surface  (Folias,  1987),  as  well  as  in  the 
interior  of  the  plate  (Folias  and  Wang  1986).  Without  going  into  the  details, 
at  z=h,  it  was  found  that 

(1) 

Og  Q 

-Qj  =  -(l+vi)  =  -  1.34,  (32) 

°rr 

suggesting,  therefore,  that  the  failure  now  is  governed  by  the  stress  Oqq'*  which 
attains  its  maximum  value  at  9  =  n/2.  Thus,  the  crack  will  begin  to  curve  into 

*  The  reader  may  notice  that  the  right  hand  side  of  equation  (29b)  differs  from  (29a)  because  it 
is  based  on  plane  stress. 
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the  matrix  until  its  direction  becomes  perpendicular  to  that  of  the  applied 
load. 


PERIODIC  ARRAY  OF  FIBERS 

The  previous  results  were  based  on  the  presence  of  one  fiber  only.  It  is 
now  desirable  to  extend  these  results  to  also  include  a  doubly  periodic  array  of 
fibers  which  are  embedded  into  a  matrix.  For  this  reason,  we  assume,  a 
periodic  arrangement  of  the  type  shown  in  fig.  9.  Following  the  same  method 
of  solution  as  that  of  Penado  and  Folias  (1989),  one  finds*  at  z=0  the  stresses 
o£jP  and  o^qP  for  vi  =  0.34,  V2  =  0.22  and  various  (G2/G1)  ratios,  shown  in  figs. 
6  and  10.  Two  observations  are  worthy  of  note.  First,  beyond  a  certain  ratio  of 
(G2/G1)  the  stress  reaches  an  asymptotic  value.  Such  trend  was  also 

found  by  Adams  (1967)  based  on  2D  considerations.  Second,  as  the  volume  of 
fibers  increases  the  stress  o^  decreases  by  as  much  as  40T  (see  fig.  11). 

Finally,  in  fig.  12  we  plot  the  stiffness  ratio  versus  the  fiber  volume  fraction 
Vf. 

Returning  next  to  the  strain  energy  release  rate,  equation  (21)  is  still  a 
good  approximation  provided  that  Go  is  replaced  by  the  following  effective 
load  stress 

f  (1)  |  1 

J  n  A  ( 

^o^ffective  periodic  '  o0)  1  '  ’  ^or  z 

[°rr  V  =  0j 

Unfortunately,  in  order  to  obtain  a  similar  expression  for  z=h,  one  needs  to 
establish  whether  the  order  of  the  singularity  strength  increases  as  adjacent 

*  The  results  are  valid  for  all  fibers  which  are  at  least  four  diameters  away  from  the  bounding 
planes  x  =  ±  w  and  y  =  ±  /  The  solution  and  the  details  are  similar  to  those  discussed  bv  Penado 

and  Folias  (1989)  except  that  one  now  has  cos(2n9),  n=  0,1,2,...,  where  the  remaining  unknown 
coefficients  are  determined  from  the  boundary  conditions  of  the  geometrical  cell 
configuration.  The  present  results,  are  based  on  n  =  0,...,  N  =  20  terms  which  provide  accurate 

results  in  the  region  |  z/h  I  <  1/2.  However,  many  more  terms  are  needed  in  oidcr  to  obtain 

accurate  results  particularly  in  the  neighborhood  of  z  =  ±  h.  We  are  presently  woiking  on  this 
and  the  results  for  this  problem,  as  well  as  for  the  problem  of  stresses  due  to  temperature 
missmatch,  will  be  reported  in  the  near  future. 
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fibers  approach  the  fiber  in  question.  In  view  of  some  previous  work  the 
author  conjectures  that  this  may  very  well  be  the  case.  Thus,  the  following 
fundamental  questions  come  to  mind.  How  close  must  adjacent  fibers  be 
before  the  order  of  the  singularity  strength  is  affected?  Does  a  certain 
separation  distance  or  a  certain  periodic  array  of  fibers  exists  which  leads  to  an 
optimal  state  of  stress?  Based  on  3D  considerations,  Penado  and  Polias  (1989) 
have  shown  that  when  fibers  are  placed  four  fiber  diameters  apart,  center  to 
center,  practically  all  fiber  interactions  have  subsided,  including  those  at  the 
free  surface  z=h.  The  author  suspects,  however,  that  when  fibers  are  placed 
two  diameters  apart,  center  to  center,  the  singularity  strength  will  be  affected. 
Naturally,  this  is  a  conjecture  that  needs  to  be  investigated. 

CONCLUSIONS 

Based  on  a  3D  analytical  solution,  we  have  shown  that  fiber /matrix 
debonding  is  most  likely  to  occur  close  to  a  free  surface.  Thus,  regions  where 
fibers  intersect  free  surfaces,  eg.  holes,  cut  outs,  edges,  cracks  etc.  are  potential 
trouble  spots.  Moreover,  the  strain  energy  release  rate  (eq.  28)  may  be  used  to 
predict  crack  initiation  in  the  center  of  the  fiber  length  (eq.  29a),  as  well  as  at 
the  free  surface  (eq.  29b).  Moreover,  fiber/matrix  debonding  at  a  free  surface 
will  occur  at  approximately  (1/10)  the  load  value  required  for  the  same  type  of 
failure  to  occur  at  the  center  of  the  fiber  length.  Such  information  on  crack 
initiation  is  particularly  important  for  the  proper  understanding  of  damage 
evolution. 

Alternatively,  the  strain  energy  release  rate  for  a  periodic  array  of  fibers 
of  the  type  shown  in  fig.  9  may,  at  z=0,  be  approximated  by  eq.  (28)  in 
conjunction  with  eq.  (33).  A  similar  expression  applicable  to  the 
neighborhood  of  the  free  surface  requires  that  one  must  first  establish 
whether  the  strength  of  the  singularity  is  indeed  affected  as  the  fiber  volume 
increases.  For  Vf  <  0.05,  however,  it  has  been  shown*  that  no  such  interaction 
effects  are  present. 

As  a  final  remark,  we  note  that  if  the  bond  at  the  interface  does  not  fail 
the  analysis  shows  that  there  exists  a  stress  magnification  factor  in  the  resin 
which  attains  a  maximum  between  the  fibers.  This  maximum  stress 


*  See  Panado  and  Folias  (1989). 
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magnification  occurs  along  the  line  0  =  0°  and  at  a  distance  r  =  1.2a  from  the 

center  of  the  fiber’* . 
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FIGURE  CAPTIONS 


Fig.  1.  Geometrical  and  loading  configuration. 

Fig.  2.  Definition  of  local  coordinates. 

Fig.  3.  Singularity  strength  for  isotropic  fiber  and  isotropic  matrix 

versus  G2/G1. 

Fig.  4.  Stress  o^at  r=a,  0=0  and  for  vi  =  0.34,  \’2  =  0.22  and  (G2/G1)  = 

16.67,  accross  the  thickness. 

Fig.  5  Stress  Gqq  at  r=a,  0=0  and  for  vi  =  0.34,  V2  =  0.22  and  (G2/G1)  = 

16.67,  across  the  thickness. 

Fig.  6  Stress  o£^at  r=a,  0=0  and  for  vi  =  0.34,  V2  =  0.22,  versus  the  ratio 

(G2/G1). 

Fig.  7  Fiber/matrix  interface  crack  under  transverse  loading. 

Fig.  8  Strain  energy  release  rate  for  plane  stress  and  plane  strain 

conditions  for  v]  =  0.34,  V2  =  0.22  and  (G2/G1)  =  16.67. 

Fig.  9  Periodic  array  of  fibers  of  length  2h,  embedded  into  a  matrix. 

Fig.  10  Stress  o^at  r=a,  0=0  and  for  v]  =  0.34,  \'2  =  0.22,  versus  the  ratio 

(G2/G1). 

Fig.  11  Stress  g^q  at  r=a,  versus  Vf,  for  (G2/G1)  =  16.67  v]  =  0.34,  V2  = 

0.22. 

Fig.  12.  Stiffness  ratio  versus  Vf  for  (G2/G1)  =  16.67,  vi  =  0.34  and  V2  = 

0.22. 
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Fig.  3. 


Singularity  strength  for  isotropic  fiber  and  isotropic  matrix 
versus  G2/G1. 


16.67.  across  the  thickness. 


a/h-0.05  r=Ji  i=Q  v  1-0.34  v2=0.22 


length  ratio  of  0.05  versus  Gf/Gm. 


STRAIN  ENERGY  AND  ITS  RELEASE  RATE 


Fig.  8  Strain  energy  release  rate  for  plane  stress  and  plane  strain 

conditions  for  vi  =  0  M.  v*»  =  n  ro*  -  tr.  m 


a/hx0-05  r= a  2=0  vl  =Q.34  v 2  =0.22 


Fig.  10  Stress  Og^at  r=a,  0=0  and  for  vi  =  0.34,  V2  =  0.22,  versus  the  ratio 
(G2/G1). 


Stress  oJJP  at  r=a,  versus 
032. 


Fig.  11 
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